General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



WASHINGTON UNIVERSITY 
DEPARTMENT OF PHYSICS 
LABORATORY FOR ULTRASONICS 
St. Louis, Missouri 63130 

(NASA"'C r? ' , -1734 16) NCN-rSS TFUCTIVE EVALUATION N84-22969 

OF COMPOSITE MATERIALS USING ULTRASOUND 
Semiannual Progress Report, 15 Sep. 1983 - 
15 Mar. 1984 (Washington Univ.) 32 p Unclas 

HC AOVtfF A0 1 CSCL 14D G3/38 00724 

"Non-Do struct! ve Evaluation of Composite Materials Using Ultrasound" 


Semiannual Progress Report: September 15, 1983 - March 15, 1984 


NASA Grant Number: NSG^lfiOl 


Principal Investigator: 

Dr. James G. Miller 
Professor of Physics 


The NASA Technical Officer for this grant is: 

Dr. Joseph S. Heyman 

NASA Langley Research Center 

Hampton, Virginia 



I 

1 

1 

a 


- 2 - 




I. INTRODUCTION 

This report summarizes our continuing investigation of the non- 
destructive evaluation of advanced composite- laminates. Oar approach 
makes use of indices derived from the measurement of fundamental acous- 
tic parameters in order to quantitatively estimate the local material 
properties of the laminate. The following sections describe ongoing 
studies of phase insensitive attenuation measurements, and discuss 
several phenomena which may influence the previously reported technique 
of polar backscatter. A simple and effective programmable gate circuit 
designed for use in estimating attenuation from backscatter is described 
in Section TV. 

II. Transmission Studies 

Previous reports from this laboratory have described experimental 
systems for the measurement of the attenuation in advanced composite 
materials by transmission of ultrasound. The procedure for phase- 
insensitive measurement of the attenuation in transmission has now been 
implemented on a LSI 11/23 computer system described in previous 
reports. The entire measurement process can now be selected and moni- 
tored from a computer terminal in the laboratory. As the data irom each 
site is recorded, the slope and intercept of the attenuation versus fre- 
quency curve is calculated and displayed for the operator's information. 
The use of the 11/23 computer system greatly speeds the acquisition and 
analysis of transmission data. 

Current studies are under way to demonstrate that the transmission 
techniques which we described previously are reproducible, 
investigato?, -independent, and quantitative. In addition, we are engaged 
in a further study of the large body of data acquired from transmission 
measurements performed on glass-epoxy composites in order to clarify the 
significance of the intercept of the attenuation versus frequency curve, 
especially in regard to (more conventional) single frequency imaging as 
compared to slope imaging. 
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III. EFFECTS OF MODE CONVERSION 

The potentially significant tool for non-des tractive evaluation 
provided by the polar backscatter technique requires a re-examination of 
the phenomena of mode conversion of sound waves at liquid/solid inter- 
faces. As we reported previously, insonifying a composite laminate at 
nonperpendicular incidence allows the selective interrogation of inter- 
nal structures, relatively independent of surface reflection effects. 
The orientation of a polar backscatter coordinate system is shown in 
Figure 1, in which the polar angle 0, and the azimuthal angle 6 are 
def ined. 

Any configuration in which an initially longitudinal wave in a 
liquid is incident on a solid at a nonperpendicular angle leads to a 
redistribution of power between the reflected longitudinal mode and the 
two transmitted modes, longitudinal and transverse. An analysis of this 
redistribution is required for proper interpretation of the energy back- 
scattered from internal structures. 

The reflection and transmission coefficients for liquid-solid mode 
conversion of plane waves in homogeneous media are available in the 
1 iterature. 1 *2 We have verified that the results obtained by these 
authors are equivalent after appropriate algebraic and trigonometric 
manipulation. We denote the pertinent angles by 0^ for the reflected 
mode, 0 2 for the transmitted longitudinal mode, and y for the transmit- 
ted transverse mode. Theso angles are related by Snell's law, as shown 
in Figure 2, The power in each mode can be represented as follows. 

2 2 "Y 

Z, co $> 2y+Z. sin 2y -Z 

R = — 5 

Zj cos 2 * 4 2y+Z t sin 2 2y +Z ± 

4 Z^ Zjcos 2 2y 

^ Z^ cos 2 2y + Z t sin 2 2y + Z^ 2 
4 Z 1 Z t sin 2 2y 

1 [z x cos 2 2y + Z t sin 2 2y + Z ± J 2 

where R is the reflection coefficient for longitudinal waves, T^ and T t 
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Figure 1. Geometry for measuring polar backscatter. The polar 
angle (0) and the azimuthal angle (0) are measured 
as shown. In this figure only one fiber orientation 
is shown; the azimuthal angle is always measured 
relative to the same fiber orientation. 
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Figure 2 — Angles of mode converted waves according 
to Snell's Law. Epoxy properties used with 0^ - 25° 


are the transmission coefficients for longitudinal and transverse waves, 
respectively, and Z represents an acoustic impodance defined by 


z i cose i 

When 9 exceed. e critical , the longitudinal mode become, imaginary, 
and is sometimes described as an "inhomogeneous wave". Letting 

cose^ s i|oos 02 l and complex conjugating appropriately, the power in the 
transverse and reflected mode are easily reformulated from the abovo 
expressions. 

Using published values for velocities in a typical epoxy, 3 we cal- 
culated the distribution of power for plane wave mode conversion (Figure 
3). The exchange of energy from the transmitted longitudinal mode into 
the transmitted transverse mode is apparent as 0 increases from zero 
degrees. Near the longitudinal critical angle of 32.4 degree, the 
reflected mode dominates. For 0 > 33 degrees, the shear mode dominates. 

Clearly this formulation is limited in several regards. First, no 

account has been taken of the energy converted into Rayleigh surface 

waves, a mode which is certainly important when 0 > 0 i . 

critical 

Secondly, the formulation does not deal explicitly with frequency depen- 
dent phenomena. (Ngoc and Mayer^ have included frequency depenaent 
attenuation effects for monochromatic reflection and transmission.) 
Finally, physical transducers do not produce "pure" plane waves, even in 
the focal zone. A useful representation is a superposition of an angu- 
lar distribution of plane waves. 2»4 This superposition accounts for the 
finite beam width, and nearly plane wave behavior in the focal zone. 
Some of the implications of these considerations will be discussed later 
in this report. 

In order to relate the considerations above to practical ultrasonic 
measurements, we used a 1.2 mm diameter piezoelectric receiver to deter*^ 
mine the beam profile of a 5 MHz focused transducer before and after a 
plastic specimen was inserted in the beam. The receiving transducer was 
scanned in a line perpendicular to the transmitter under control of the 
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Figure 3 — Distribution of power for mode conversion at water/epoxy interface. 
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PDP 11/23 computer system described in previous progress reports. A 
selected portion of the received signal vas gated into a spectrum 
analyzer. The 6 MHz component was read by the computer ana stored on 
disc for later processing, A 5 mm thick slab of Luoite was then intro- 
duced between the transducers, which were mounted to allow rotation of 
the slab to achieve various incxdont angles. The thickness of the slab 
was chosen based on handbook5 values of the velocities in Luoitot, Sim- 
ple calculations showed that this thickness would result in a sufficient 
time separation of any transverse waves from the faster longitudinal 
waves. 

Figure 4 comperes the beam profile measured with a water-only con- 
figuration to the beam profile obtained after transmission through the 
plastic at normal incidence. The horizontal scale on this and the sub- 
sequent four figures represent the lateral dimensions of the scan line. 
Although the absolute origin of the horizontal axis varies in some fig- 
ures, the scale size is unchanged. The backscatter was calibrated in dB 
relative to the central peak of the beam. The water-only profile was 
calibrated relative to itu peak. All other profiles were calibrated in 
terms of the central peak of the normal incidence ( 0 * 0 degrees) 
transmitted profile, which was 58 dB below that of the water-only pro- 
file. For convenience, the 58 dB difference is omitted in Figure 4. We 
note that after the beam was transmitted through the plastic plate at 
normal incidence, the relative shape of the beam was essentially 
unchanged from that exhibited by transmission through a water-only patn, 
except that the sidelobes are about 5 dB higher(with respect to the mam 
lobe) after insertion of the Lucite. 

Aft^r a rotation of the sample to an incident angle of 50 degrees, 
the critical angle of 34 degrees is definitely exceeded; no signal is 
observable at the arrival time of any longitudinal modes. Thus, the 
beam profile of the transverse mode should contain almost all the energy 
transmitted into the sample. In Figure 5 the pure L-wave profile shewn 
in Figure 4 is compared with the pure T-wave observed at 50 degrees. To 
facilitate comparison, the profiles have been aligned horizontally, 
since the shear wave experiences refraction. The profiles are similar, 








with the sidtlobes on the she «r wave located In the sane relative posi- 
tion ai in the water-only cate. The difference in magnitude of these 
nodes nay be related to the energy converted into surface wave modes, as 
well as differences in the attenuation of each node. 

The exchange of energy between the transmitted longitudinal and 
transverse nodes as the incident angle is varied is depioted in Figures 
6 , 7, and 8. The angles of incidence are 20, 3C> and 40 degrees, 
respectively. The horizontal scale on each figure is self consistent, 
so that the actual separation of the beans can be measured. We note 
that near the critical angle (Figure 7) and even just beyond the oritioal 
angle (Figure 8), there is appreciable energy arriving at the appropriate 
time for longitudinal waves. This ''smearing*' of the longitudinal peak, 
as well as the existence of energy in longitudinal waves beyond the 
critical angle support the view of a finite ultrasonic bean as a super- 
position of an angular dlstribti^ ion of plane waves. Thus for any angle 
of incidence there are some plane waves, albeit small in magnitude, 
which do not exceed the critical angle. When the majority of the plane 
waves were near the critical angle, where small variations in incident 
angle result in dramatic variation in transmitted angle, these waves 
experience refraction, which results in an angular distribution of 
energy. 

Conversely, this angular distribution interpretation suggests that 
there is always some amount of energy incident at or beyond any of the 
critical angles of the material. The next section of the report 
discusses the implications of one such angle, the Rayleigh critical 
angle. 

TV. Bounded Beam Surface Effects 

The scattering of a bounded acoustic wave from a smooth liquid- 
solid interface has been studied in the past few years by several inves- 
tigators. One most interesting effect described is a relative maximum 
in scattering at the Rayleigh angle. This result was reported by de 
Billy, Adler, and Qnentin<> for water-stainless steel and water-copper 
interfaces. These same investigators later observed a relative maximum 


6 MHz Beam Profile, 20 Degree Incidence 
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Figure 6 - Comparison of mode converted beam profiles, demonstrating the 
distribution of power and horizontal displacement at 20° degree incidence 







6 MHz Beam Profile, 40 Degree 
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at bouh the forward and backward scattering Rayleigh angles, for any 
incident angle. 7 Several theoretical methods have been investigated to 
attempt to explain the experimental reports. Although elegant formal 
presentations have been ueveloped8,9,10 these are not as computationally 
straightforward as a numerical algorithm developed by Ngoc and Mayer. 4 
All of the proposed theories appear to agree in the prediction that the 
relative naximum in backscatter at the Rayleigh angle must be related to 
the bounded aspect of the acoustic beam as well as to the resonant Ray- 
leigh surface propagation mode. 

The possible implications of this effect for our studies of fatigue 
and impact damage in composite laminates prompted us to undertake a 
short study of the effect. A planar stainless steel plate was insoni- 
fied by a 5 MHz center frequency transducer using a broad band width 
pulse, as opposed to the single frequency studies reported in the 
literature. The incident angle was swept through 60 degrees, in one 
degree increments, under computer control, with the entire backscattered 
signal gated into a computer controlled spectrum analyzer. Several fre- 
quency components were recorded onto disc at each one degree increment. 
Results were calibrated with respect to the normally incident reflec- 
tion. Figures 9 and 10 show the 5 MHz component of direct backscatter 
between 10 and 60 degrees. The vertical scale of Figure 9 is linear, 
that of Figure 10 is logarithmic. We note the distinct peak centered at 
approximately 32 degrees. The predicted Rayleigh angle for stainless 
steel, based on a velocity of 2891 m/sec, is approximately 31 degrees. 
The magnitude of the relative peak is approximately 27 dB below that of 
normal incidence. This magnitude, as well as the angular spread of the 
peak, is thought to be related to the size of the transducer ana the 
spatial width of the ultrasonic beam. Ngoc and Mayer predict that a 
larger diameter transducer and tighter focusing result in narrowing the 
width of the peak, with transducer width paralleling variations in the 
relative magnitude, and tightness of focusing inversely related to the 
relative magnitude. 

Thftse considerations provide an approach to measuring the Rayleigh 
angle and hence the Rayleigh velocity. Of more direct interest to our 
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5MHz Polar Backscatter From Steel 



m 

r d 


Sh 

o 

•+j 

o 

a) 

r — H 

M-t 

<u 

K 

a 

0 

u 

0 

Oh 

£ 

O 

i — i 
0 
PQ 


0 

— 10 )— 

t 

j 

— 20 i — 

— 30 

— 40 


Logarithmic Scale 


-50 


0 
-P 

cd 
o 

J 3 -60 

o 
cd 

pq 



0 


10 TO 2 ! 0 A pO / r . 40 ' 50 

Polar Angle (Degr-'^s) 


60 
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efforts, however, it the relationship of these considerations to the 
polar backscatter measurement* that we previously reported. Polar back- 
scatter in graphite epoxy composites subjected to impact and fatigue 
damage was between 25 and 35 dB below that of a pertect (perpendicularly 
oriented) reflector. Clearly, if the insonifying angle was close to the 
Rayleigh critical angle, the bounded beam effect described above would 
have the potential to mask scattering from internal structures. 

Further investigation of this effect may clarify the origin of dif- 
ficulties in using polar backscatter on glass/epoxy specimens described 
in the previous report. For example, the Rayleigh velocity in fused 
silica glass, as determined from the Rayleigh wave velocity equation, H 
is 3400 m/sec, which leads to a Rayleigh critical angle, in water, of 
approximately 26 degrees. Therefore, if a glass-epoxy laminate was suf- 
ficiently glass-like nltrasonically, especially in terms of shear wave 
velocity, then the backscatter at the Rayleigh angle could corrupt meas- 
urements taken at a polar angle near the critical Rayleigh angle. 

We are planning to conduct additional investigations of the effect 
of broadband vs narrow band insonif ication, as well as of the effects of 
velocity dispersion. 

V. Anisotropic Mode Conversion 

Considerations developed in the previous section suggest that the 
Rayleigh critical angle (and possibly the longitudinal critical angle as 
well) should be avoided when using polar backscatter techniques. In 
isotropic media these angles can be calculated or measured straightfor- 
wardly. In contrast, in anisotropLs materials such as composite lam- 
inates the values of these critical angles are markedly dependent upon 
the polar and azimuthal angles of the insonifying transducer, as defined 
in Figure 1. 

The division of power between the modes is very important for the 
interpretation of quantitative nondestructive evaluation, especially in 
the calculation of quantitative backscatter coefficients, where the spa- 
tial energy distribution of the beam must be folded into the 
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calculation. 12 The calculations in this section will be based on experi- 
mental velocity values in unidirectional graphite-epoxy composites as 
reported by Williams. 3 The Rayleigh angle backscattering maxima dis- 
cussed in the previous section will be ignored, since the Rayleigh angle 
is estimated to be greater than SO degrees, which is well beyond the 
range of polar angles we will consider here. 

Williams and Nayeb-Hashemi measured the velocities of ultrasound in 
orthogonal directions in unidirectional graphite epoxy. As one would 
anticipate the velocities were fastest for propagation parallel to the 
fibers and slowest for propagation perpendicular. For example, the 
longitudinal velocity parallel to the fibers was 8400 m/sec and the per- 
pendicular longitudinal velocity was 2800 m/seo. In general, the 
transverse velocities depend on the polarization as well as the propaga- 
tion direction, varying between 1270 m/sec for the slowest perpendicular 
orientation ( both propagation and displacement perpendicular to fibers) 
to 2000 m/sec for the parallel orientation (propagation parallel to 
fibers, displacement perpendicular to the fibers ). 

A detailed study of the exact functional dependence of the veloci- 
ties on 6 and i in such an inhomogeneous, anisotropic medium is not con- 
sistent with the goals of the current research. However, some insight 
may be gained into the distribution of power in mode conversion by mak- 
ing a few simple approximations. First we assume that the velocities 
vary smoothly as 0 and i vary, perhaps sinusoidally. Secondly, we 
assume that the plane wave approximation holds so that the coefficients 
developed in section three can be used by replacing the constant veloci- 
ties with the appropriate function of 0 and 6. 

The velocities calculated up*’er the first assumption are depicted 
in Figures 11 and 12 . Figure 11 contains families of curves represent- 
ing the transverse velocity at fixed polar angles plotted as functions 
of azimuthal angle. The curves are shown for polar angles of 5 , 10 , 15 , 
20 , 25 , and 30 degrees. Although the azimuthal variation is signifi- 
cant, the variation with polar angle is less dramatic, since the range 
of polar angles is well away from the critical angle for shear waves. 


Transverse Velocity in Unidirectional Composite 







- 21 - 


Figure 12 depicts similar families of curves corresponding to the longi- 
tudinal velocity for polar angles of 5 degrees through 30 degrees as in 
Figure 11. As the polar angle increases, the range of azimuthal angle 
over which the longitudinal wave exists decreases resulting in the pro- 
gressively narrower constant 8 curves in Figures 12 through 16. This 
result is illustrated in Figure 13 in which the longitudinal critical 

angle is plotted as a function of d. Clearly this variation in 0 , it 

criti- 
cal will have a significant impact on the proper design of polar back- 
scatter measurements. 

We now apply the second assumption and utilize these velocities in 
the power coefficients from Section III. For simplicity we ignore the 

longitudinal mode for 0 > under the reasonable approximation 
that most of the power will be in a transmitted transverse mode for 
those angles. Figures 14, 15, and 16 display families of curves for 
polar angles of 5 degrees through 30 degrees in 5 degree increments used 
in Figures 11 and 12. Each curve represents the fractional power in 
that mode, for the specified polar angle , as a function of azimuthal 
angle. Figure 14 depicts the power in the reflected longitudinal mode, 
which exhibits considerable variation with d, especially near the criti- 
cal longitudinal angle for that d. Figure 15 plots the power in the 
transmitted longitudinal mode. Figure 16 displays the power in the 
transmitted transverse mode. Both transmitted modes show substantial 
variation near the critical angle. 

The results of these elementary calculations have several implica- 
tions. First the amount of power reflected from the surface depends 
significantly on the polar and azimuthal angles especially near the 
local critical angle. This signal can contribute significantly to the 
total polar backscatter. This combined with the implications of the 
Rayleigh angle effect described in section IV suggests that the polar 
backscatter technique in composites is not totally free of surface 
effects. We note, however, that with care, these effects can be minim- 
ized. Further consideration of improvements to the polar backscatter 
technique in order to reduce surface effects is under way. A second 
implication is the dramatic variation in the longitudinal critical 



















angle. For 6 > ® cr ^ioal' wo$t power is in the transverse mode, 
with some power in surface waves and more in the longitudinal mode. The 
two regimes defined by 6 < 0 Q and 0 > 0 C can be achieved at the same 
polar angle merely by rotating the azimuthal angle. Thus the calcula- 
tion of the bachscatter coefficients from the measured backscatter 
transfer function requires special attention under these conditions. In 
addition, the scattering properties exhibited by zones of damage differ 
depending upon the interrogating sonic mode. We plan further investiga- 
tions of the effects of mode conversion in anisotropic media and the 
scattering properties of the different modes. 

VI, Programmable Gate Circuit 

Work in this and other laboratories has demonstrated the potential 
value of quantitative measurements of the slope of the attenuation for 
the nondestructive evaluation of composite materials. The transmission 
techniques used to date may not be feasible in many circumstances, lead- 
ing to a demand for techniques of estimating the attenuation from 
reflected ultrasound, A number of such techniques have been pro- 
posed. 13, 14, 15, 16 

One such technique involves measuring the shift of the centroid of 
the frequency spectrum of the insonifying beam as a function of deptn in 
the attenuating medium. This approach requires a knowledge of the fre- 
quency spectrum from a number of depths, obtained either by digitization 
of the time domain signal and subsequent processing, or by gating 
appropriate segments of the signal into an analog spectrum analyzer, 
from which the desired spectra can be recorded directly. 

Figure 17 is a sketch of a typical measurement. The symbol W 
represents the width or duration of the gated segment ana D is the 
overall time delay to the beginning of the gate. The width of the gate 
is selected and the delay moved through the sample so that the spectrum 
from each depth can be recorde We intend to use this system to esti- 
mate the slope of the attenuation in a very thick and highly attenuating 
sample of graphite/epoxy composite cut from a tag-end of the filament 
wound casing being tested for the Shuttle solid rocket motor. 
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The circuit to achieve this programmable gate it depicted by the 
block diagram of Figure 18* The circuit was constructed on a Hewlett- 
Packard Multiprogrammer Breadboard Output Card(model 69380A), This 
breadboard card stores 12 bits of data from the computer using built-in 
data latches. The circuit can be visualized as three stages built 
around a timing source, The timing source for this application is a 
20,000 MHz K1114AM crystal oscillator. The three stages consist of a 
time delay counter, a gate width(or duration) counter, and an analog 
gate output. 

The timing of the delay counter is derived from the 20 MHz clock by 
using a 74LS169 synchronous 4-Bit up/down counter no a divider, with the 
rate determined by a 4“bit DIP switch mounted on the board. The avail- 
able delay clock rates are then between 1.3 MHz and 10 MHz, For most 
applications in this laboratory, 10 MHz is the preferred rate. This 
clock rate drives a 12 bit delay counter composed of three 74LS169 chips 
set up in cascading 4-bit stages. The counter is enabled and commences 
counting down by an external trigger, which triggers a 74123 one-shot 
into the J-port of a 74H73 high speed JT-K flip flop. When the delay- 
counter finishes, the carry out pulse triggers two 74123 one-shots, one 
of which reloads the delay data from the latches, while the other tog- 
gles the K input of the delay counter flip-flop, terminating counting. 
This termination pulse is also directed to the J input of a second 
flip-flop, which initiates the gate width counter and serves as the 
source of a TIL gate of the appropriate delay and duration. 

The gate duration is determined by 8 bits run at 20MHz. The 8-bit 
information can be entered by hand via DIP switches on the board, or 
programmed via a dfgital output card. The gate width counter consists 
of two 74LS169 chips in cascade, much as the delay counter was con- 
structed. The carry output, upon completion of the countdown, triggers 
a one-shot which reloads the width data and toggles the K input of the 
flip-flop, disabling the counter and terminating the TTL gate. 

The TTL gate requires some analog shaping and buffering. The Q 
output of the J-K flip-flop is run through a LM318 Op-Amp set up as an 
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inverting discriminator* driving anything below 1.5 volts("on") to 15 
volts and anything above l,5volts ("off") to -15 volts thus eliminating 
"ringing" on the TTL signal. The resulting signal is clipped to approx- 
imately 5 volts using a Zener diode and the resulting gate is buffered 
through a 2N2222 transistor to supply the necessary drive levels. 

The oircuit so described provides a nicely shaped gate between 0.6 
and 12.4 microseconds in duration, with a maximum delay (at 10 MHz) of 
409 microseconds. The software required to drive the gate is very sim- 
ple and allows quick and accurate placement of the gate under control of 
our LSI 11/23 or HP9825 computers. 

VII. Conclusions 

The main thrust of this report has been the discussion of some sur- 
face effects which may influence polar backscatter measurements, and the 
implications of mode conversion in calculating a quantitative back- 
scatter coefficient. Further study of the details of mode conversion in 
inhomogeneous anisotropic media is needed to enhance our continuing 
study of the backscatter coefficient in composite-laminates. The bound- 
edness of the insonifying beam may prove troublesome for polar back- 
scatter in some media, especially near the Rayleigh angle. With some 
knowledge of the material under investigation, the troublesome effects 
of the bounded beam can be minimized by selecting polar angles away from 
the critical Rayleigh angle. 
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